Hopping in the Glass Configuration Space: Subaging and Generalized Scaling Laws 
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Aging dynamics in glassy systems is investigated by considering the hopping motion in a 
rugged energy landscape whose deep minima are characterized by an exponential density of 
states p(E) = Tg 1 exp(E/T g ), — oo <E<0. In particular we explore the behavior of a generic 
two-time correlation function TL(t w +t, t w ) below the glass transition temperature T g when both 
the observation time t and the waiting time t w become large. We show the occurrence of ordinary 
scaling behavior, IL(t w +t,t w ) ~ Fi(t where fxi = l (normal aging) or /ii < 1 (subaging), 
and the possible simultaneous occurrence of generalized scaling behavior, t w [l — Tl(tw +t, tw)] ~ 
Ft2(t/t w 2 ) with fi2 <Mi (subaging). Which situation occurs depends on the form of the effective 
transition rates between the low lying states. Employing a "Partial Equilibrium Concept" , the 
exponents jiti,2 and the asymptotic form of the scaling functions are obtained both by simple 
scaling arguments and by analytical calculations. The predicted scaling properties compare well 
with Monte-Carlo simulations in dimensions d = 1 — 1000 and it is argued that a mean-field 
type treatment of the hopping motion fails to describe the aging dynamics in any dimension. 
Implications for more general situations involving different forms of transition rates and the 
occurrence of many scaling regimes in the t-t w plane are pointed out. 

PACS numbers: 02.50.-r, 75.10.Nr, 05.20.-y 



I. INTRODUCTION 

The history of glass formation strongly affects the re- 
laxation dynamics of glassy materials (l], |^]. This dy- 
namics is found to become slower with the "age" of the 
system, that means with the time t w expired since the 
material was brought into the glassy state. Such aging 
phenomena have been identified in many systems and 
various dynamical probes. Prominent examples are shear 
stress relaxations in structural glasses ||, thermorema- 
ncnt magnetizations or a.c. susceptibility in spin glasses 
0, p. Similar effects have been observed on the di- 
electric constant of dipolar glasses || 0], of structural 
glasses ||, [|, and on the structure factor of Lennard- 
Jones systems |l0| . More recent experiments in colloidal 
gels [[[T or other 'soft glassy materials' have been re- 
ported[p|, pi, II @ , and even electronic relaxations 
in Anderson insulators were found to exhibit aging effects 
|l7f . Aging is also expected for pinned systems (pinned 
domain walls pinned vortex lines fTil ), polymer 

melts ||(J and granular materials (see e.g. |pl|, p2[). 

From a theoretical point of view, several pictures have 
been proposed p3|. The simplest one is based on domain 
coarsening ideas |24|, and is probably well suited to de- 
scribe aging in, say, disordered ferromagnets where a well 
defined order wants to establish across the system. How- 
ever, in spin glasses and even more evidently in glasses, 
the idea of some long range order which progressively in- 
vades the system is far from trivial. Mean field models for 
spin glasses, which are formally equivalent to the Mode 
Coupling Theory of glasses p3[], do indeed lead to aging 



phenomena below the glass transition. In this case, aging 
is of geometric origin p5fl : As time grows, the system pro- 
gressively exhausts the possibilities of lowering its energy, 
and finds itself around saddle points from which it is more 
and more difficult to escape. Thermal activation is irrel- 
evant in these models. Although this picture might be 
justified for supercooled liquids above the Mode-Coupling 
temperature [ p6| , it certainly breaks down at lower tem- 
perature, where activated events become dominant. In 
this regime, one expects that a coarse-grained dynamical 
model of thermally activated hops between metastable 
states is a proper description of the dynamics. In fact, 
recent MD simulations support this view [g7| p9[ . 
Landscape models have been widely discussed in the past 
p(j| . |3l| , but their relevance for aging phenomena was rec- 
ognized later 



p5| . These later developments 
were recently extended to treat rheological phenomena 

s 

The 'trap' models studied up to now lead to correlation 
or response functions that depend on the ratio t/t w of the 
observation time t to the waiting time t w (full aging), or, 
for long-ran ge c orrelated energy landscapes, on the ratio 
logi/logiu, |37|| . However, many experimental systems 
reveal subaging behavior, that is, the relevant variable is 
t/t^ with fi<l. Furthermore, it is possible that there ex- 
ist, for given waiting time t w , various scaling regimes in 
time t, which are governed by different relaxation times 
oc t^ E , s — 1,2,.... The occurrence of different scaling 
regions has recently been conjectured on the basis of an- 
alytical results for mean-field spin glass models |2^, p8[ . 
So far, however, it was not possible to illustrate intu- 
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FIG. 1: Sketch of a rugged energy landscape with various 
metastable minima. Within a coarse-grained description, the 
slow dynamics of the system may be attributed to effective 
"super-transitions" between the deepest minima belonging to 
the shaded area. 



itively these multiple time regimes by exact calculations 
on simpler models (see however |39| for an interesting 
discussion of these multiple time scales). In this paper 
we will discuss a model that allows us to demonstrate for 
the first time explicitly the possible occurrence of sub- 
aging behavior and multiple time scaling in a hopping 
model, where a point jumps among the deep (free) en- 
ergy minima of a complex configuration space. Some 
of the results discussed in this paper already appeared in 
a Letter El. 



II. HOPPING IN A RANDOM ENERGY 
LANDSCAPE 

A. Model 

Slow dynamics in glasses is often attributed to a ther- 
mally activated motion of a point (henceforth denoted 
as "particle") that jumps among metastable states in 
a rugged configuration space (see Fig. p. In a coarse 
grained description, only transitions between the deepest 
(free) energy minima Ei will govern the dynamics at long 
times. According to extreme value statistics one may ex- 
pect the distribution ip(E) of these deep minima to be 
exponential (which is the behavior of the tail of a Gumbcl 
distribution, see e.g. [flO)). Indeed, mean- field theories of 
spin glasses Gj] and recent results from molecular dy- 
namics simulations |42| suggest this to be the case. 

For simplicity, we consider the metastable states with 
lowest energies to be arranged on a hypercubic lattice in 
d dimensions and refer to them as "sites". The lattice 
here resembles an average finite connectivity of the mu- 
tually accessible states. To each lattice site i is assigned 
a random energy Ei drawn from a distribution 



il>(E) 



I ( E 

exp 



-oo < E < 



(1) 



As discussed in a moment, T g corresponds to a "glass 
transition temperature" and we thus define 



9 ee T/T g 



(2) 



as the rescaled temperature. Also, energies are specified 
in units of T g . The particle can jump from one site i to 
any of the 2d nearest neighbor sites j with a hopping rate 



v exp 



[aEj - (1 - a)Ei] 



(3) 



where the "attempt frequency" is = l sets our time unit, 
and the parameter a specifies how the energies of the ini- 
tial and target site contribute to the saddle point energy 
being surmounted during a jump. In order for the w% j to 
obey detailed balance, a can assume any real value, but 
in a sense of a weighting of the initial and target energy 
we restrict a to the range 



< a< 1 



(4) 



Note that the case a = defines a trap model, where the 
jump rates depend on the initial energy only, while the 
case a = 1/2 defines a "force model", where the jump 
rates are determined by the energy difference between 
the two sites. 

Due to the existence of very deep traps, the system ex- 
hibits a "dynamical phase transition" at the glass transi- 
tion temperature 6 = 1: in the high-temperature regime 
9 > 1 the Boltzmann distribution 



(5) 



is normalizable, while in the low-temperature an equi- 
librium state does not exist. In this latter situation the 
system is exploring deeper and deeper traps as time pro- 
ceeds and the overall dynamics ages. 

For the trap model (a = 0) the dynamics is fully char- 
acterized by the trapping times 



exp(E/0) , 



with distribution 



P (t) = 9t 



1 < r < oo 



(6) 



(7) 



The absence of an equilibrium state for 9 < 1 is reflected 
in the fact that the mean trapping time (r) becomes in- 
finite. For a > it is convenient to operate with the r 
defined in (||) as well, although these r no longer have 
the meaning of a trapping time. The hopping rate (||) 
can then be written as 



u *,3 



(8) 



At first glance the model defined here seems to be sim- 
ilar to the model considered earlier in j35|. However, 
there is an important difference. In the model studied in 
p5| the disorder is of "annealed" type, i. e. the energy 
of each site is drawn anew from tp(E) after each jump. 
This can be viewed as a mean-field treatment and sim- 
plifies the problem to a great extent. As will be shown 
in sec. HI, in the annealed case one can map any param- 
eters (9, a > 0) to parameters (6', a — 0), that means 
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FIG. 2: Aging function U(t + t w ,t w ) for (d,0,a) = 
(10, 1/4, 3/8) and different waiting times t w . The symbols re- 
fer to waiting times t w = 2 x 10 8 (+), 5 x 10 s (o), 2 x 10 9 (♦), 



5xl0 M (□), lxlO 10 (■),4xl0 1 
and 1 x 10 12 (a), respectively. 



(o), lx 10 11 (•), 3x 10 11 (a), 



the a parameter is irrelevant. In the present model by 
contrast the energy disorder is "quenched" and such a 
mapping is not possible. We will show that this leads to 
the occurrence of much richer aging dynamics, including 
subaging behavior and generalized scaling forms. 



B. Aging function 

For studying aging properties we consider a quench 
from 9 — oo to 9 < 1 at time 0, then wait for some time 
t w 3> 1 and ask for the behavior of correlation functions 
followed during an observation time t after the waiting 
period. Since the physical observables are functions of 
the coordinates of the configuration space, they will es- 
sentially de-correlate after the system has undergone a 
single transition from one deep minimum to another one. 
Hence, a "generic correlation function" in the model is to 
consider the probability that the particle has not jumped 
between t w and t w + 1. We denote this probability, after 
performing the disorder average, by U(t w +t,t w ). 

Let us note that Tl(t w +t,t w ) can also be regarded Q 
as the spin-spin-correlation function of the generalized 
Sherrington-Kirkpatrick model with p-spin interac- 
tions in the limit p — > oo, which maps onto the random 
energy model Ej| . It can also be viewed as an incoherent 
intermediate scattering function at large wave numbers 
in diffusion processes M. 

A typical result for Tl(t + t w , t w ) from continuous-time 
Monte Carlo simulations is shown in Fig. for parame- 
ters (d,9 7 a) = (10,1/4,3/8) and different t w (for details 
regarding the simulation see appendix |A|). Indeed, we 
can identify a pronounced aging phenomenon: The de- 



cay of Tl(t w +t, t w ) with time t becomes slower and slower 
with increasing waiting time t w . The linear behavior of 
the curves at long times in the double-logarithmic plot 
indicates a common asymptotic power-law decay, and 
suggests that the various curves may be collapsed onto 
a common master after rescaling the observation time t 
by a proper function of the waiting time t w . The under- 
standing of such scaling properties of U(t + t w , t w ) when 
both t and t w become large (t, t w 3> 1) will be the central 
issue in the following. 



III. ANNEALED SITUATION 

In the annealed case the model corresponds to a 
continuous-time random walk (CTRW) irrespective of 
the value of a, since all times Tj are drawn anew after 
each transition. The whole process can be characterized 
by a distribution p e g(T c g) of effective trapping times r e g 
(residence times between two transitions). 

In order to calculate p c g{T e g), let us consider the par- 
ticle to be at site i = 0. Then 



1 



Toff 



2,/ 



2d 



(9) 



and the distribution is given by 



Pefr(feff) 
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l/(l-o) 
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(10) 



The S function contributes if T e g Ej=i T f — 1> i- e - there 
is no restriction for r e g > l/2d. Hence, for r e g > l/2d: 



Pcff(Tcff) 



C e ff6'T e g~ 



e' 



(i 



(ii) 



where C e g = (ELi t > 



3 = 1 '3 



>, and (...) denotes an aver- 
age over 2d uncorrelated random numbers Tj distributed 
according to (^). The CTRW theory applies with a 
rescaled temperature 9' ', and the aging properties of this 
model have been worked out in pa ]. 

Note in particular that for 6^~> 1 no aging occurs, 
even if the distribution (^) is not normalizable for 9 < 1. 
This, however, is nothing to worry about, since in the an- 
nealed model the energies change after each hop and the 
distribution (|^) does not correspond to an equilibrium 
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distribution. The dynamical phase transition in the an- 
nealed model is defined by the diverging mean waiting 
time of the distribution p c s(t c s) that occurs for 6' < 1. 
The absence of aging for < 1 — 8 < a, that is when 
9' > 1 but < 9 < 1, points to the fact that the annealed 
model may not provide a valid mean field description of 
the quenched model in any dimension d. 

In Fig. H the aging function H(t w + t,t w ) is shown 
for the annealed model as resulting from simulations 
for two different parameters sets, (d,9,a) — (10,1/3,0) 
and (d,9,a) = (10,1/4,1/4) that both refer to the 
same 9' — 1/3. As expected, the aging functions for 
the two parameter sets are indistinguishable. More- 
over, as predicted in Q, U(t w + t, t w ) scales with t/t w , 



Tl(t w +t,t w ) = F{t/t w ), where F(u) ~ 1- 



for u <C 1 



(see Fig. |^) and F(u) ~ u~ x (see Fig. |3ja). We will see 
in the following that this simple behavior does not hold 
true any longer in the quenched case. Nevertheless, since 
T e g is the inverse hopping rate on a site, the formulae 
(plO) will still be useful in the following. 



IV. PARTIAL EQUILIBRIUM CONCEPT (PEC) 

Although full equilibrium cannot be reached in a sys- 
tem of infinite size, there should be some equilibration 
on a local scale that corresponds to the region of con- 
figuration space being explored by the system after the 
quench. This is the idea of the PEC. 

In the present model we can translate this idea into a 
precise though approximate procedure for describing the 
aging process. After the waiting time t w , the particle has 
visited S = S(t w ) distinct sites and we will assume that 
on these sites the system has equilibrated. Accordingly, 
the probability pj to find the particle on the site j of the 
set of visited sites is 



Pj (S) 



-Biie 



fc=i 



o-E k /e 



Tk 



(12) 



The probability for the particle to remain on site j for 
a time t is exp(-t£„ 3 Wj, nj ) = exp(-tr" _1 r«.), 
where the sum over rij runs over all nearest neighbor 
sites of j. Since the distinct visited sites form a Brow- 
nian path in the lattice, which, on a local scale, has a 
one-dimensional topology, we will consider exactly two 
of the rij sites to belong to the path. The remaining 
n — 2 sites are considered to be never visited until t w . 
It will turn out (see sec. VII) that this assumption is 



not very important and that the aging properties would 
be mainly the same if the Brownian path had a com- 
pact structure (except for the asymptotics discussed in 
According to the PEC, the aging function 
+ t,t w ) is now approximated by taking the average 



sec 



of exp(— t t" 1 Y^n- T rij ) over au visited sites j with the 



weights pj , 
IL(t w + 1, t v 



Egr^.exp^r;- 1 ^.^ 



(13) 



Here (...) denotes an average over (2d— 1) S(t w ) uncorre- 
cted random numbers Tj that are distributed according 
to We will see that U(t w + t,t w ) and U(t w + t,t w ) 
exhibit the same scaling properties, which in view of the 
generic character of H(t w + t,t w ) (cf. Sec. [IB ) are of 
primary interest for us here. In quantitative terms both 
functions turn out to be different. Since Tl(t w + t,t w ) 
depends on t w only through S = S(t w ), we will, in the 
following, also denote this function by Il(t, S). 

As shown in appendix |b| the number of distinct visited 
sites S(t w ) grows with increasing waiting time t w as 



S(t w ) ~ 



where 



= ) d+(2-d)6 



for d < 2 



for d > 2 



(14a) 



(14b) 



In d = 2 there are logarithmic corrections, S(t w ) 
[t w /\ogt w ] e . 



V. SIMPLE SCALING ARGUMENTS 

In this section we present simple arguments, why and 
how H(t w + t,t w ) is expected to scale. In particular 
we propose the presence of two characteristic relaxation 
times U(t w ), i — 1,2, that grow as power laws with t w , 



(15) 



where fii are exponents depending on 9 and a. 

The first characteristic relaxation time t±(t w ) can be 
obtained as follows. After the waiting time t w , the par- 
ticle has visited S — S(t w ) distinct sites and the oc- 
cupation probability of these sites is dominated by the 
one for the site with lowest energy i? m i n (in the glassy 
phase 9 < 1 the extreme values are dominant). The typ- 
ical value of -Emin can be estimated from S J_ ™ ip(E) — 
Sexp(E min ) ~ 1, i.e. 



E min (t w ) ~ - hi S(t w ) - -jhi(t w ), 



(16) 



or, in terms of the variable r max corresponding to E m i n , 
r m ax(i») = exp(-£ min /0) ~ 5(^) 1/e - t^ e . (17) 

Being at the site with energy -Emm, the particle typically 
encounters a situation as drawn in Fig. [|a: All neighbor- 
ing sites have energies close to zero and the characteristic 
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FIG. 3: Aging functions (a) H(t + t w ,t w ) and its complement (b) 1 — H(t + t w ,t w ) from simulations of the annealed model with 
parameters (d,8,a) — (10, 1/3,0) (lines) and (d,6,a) = (10, 1/4, 1/4) (points). The waiting times range from 10 9 to 10 12 . The 
scaling is so good that the various curves cannot be distinguished on the scale of the line thickness. 




neighbor 




^neighbor ^min 



FIG. 4: Jumps out of a valley with energy -E m i n after time t m : 
(a) a situation corresponding to the typical situation, where 
the energies of the neighboring valleys have values close to 
zero; (b) a rare event where at least one of the neighboring 
valleys has an energy of order -Emin. 



escape time t\ from the site with energy E mnl is the in- 
verse hopping rate given in eq. @, ti(t w ) ~ r^" — 



i.e. 



Mi 



7(1 ~ a) 
9 



(18) 



Hence, for both t and t w becoming large (t,t w 3> 1) and 
Ai = t/t!£ being fixed we expect Il(t, S) to become a 
function of Ai only, 



IL(t,S)~F 1 (A 1 ), Ai=t/C 



(19) 



The second relaxation time t2(t w ) is associated with 
rare events depicted in Fig. |Jb, which turn out to be im- 
portant at small times. The deepest state on an interval 
of length S has a typical value determined by the fact 
that it should occur with a probability 1/5. To obtain 
the scaling of the second deepest state with S we note 
that for a distribution exhibiting no peculiar long-time 
tails, the gap between the deepest and second deepest 
state remains finite when S — > oo. Accordingly, with 



probability of order 1/S ~ r~f x one of the neighbor- 
ing sites can also have an energy comparable to E m in- 
In this case the hopping rate specified in eq. (||) yields 



ti 



-a l a ^ a(l-2a)/6 
ix / max 



/'2 



7 (1 - 2a) 



(20) 



Clearly this is a relevant time growing with t w only for 
a < 1/2. Moreover, there is a further condition on the 
relevance of ti. To see this, we may consider the formal 
short-time expansion of n(i, S) from eq. ([l3|). Using the 
fact that a sum over random t's raised to some power 
with an exponent larger than 9 scales as the maximal 
term appearing in the sum (Levy statistics), we can esti- 
mate 



\m— j—l rij 



■ (1— ct)m 
m—1 ^ ■ ^max 



1 +T~ 8 T am 
1 max max 



(21) 



Here the first term in the rectangular brackets corre- 
sponds to the typical situation, while the second term 
corresponds to the rare events and is weighted by a fac- 
tor r~f x (we have neglected constant prefactors). Com- 
paring these two terms, we see that, if 9 < a, the latter 
term dominates for all m. In this case we thus expect 
the scaling behavior of the form 



S[l-fL(t t Sj\~F 2 (A 2 ), A 2 = t/t 2 (t w ) 



(22) 



for times t, t w ^> 1 and 9 < a < 1/2. 

For a < 9 to the contrary, the first term in the rect- 
angular brackets dominates for small 1 < m < 9 /a, and 
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t 2 (t w ) is not significant [note that t 2 (t w ) is smaller than 
ti(t w )]. The fact that eq. ( |2l| ) is not a regular expansion 
in the scaling variable t/t\{t w ) = £/ r max (for a > 0) in- 
dicates that Fi(Ai) does not have an analytical behavior 
for small Ai. On the other hand, the formal short-time 
expansion suggests a regular behavior F 2 {h. 2 ) ~ A2. An 
exact treatment of the PEC formula (O) presented in 
however, shows that this should be true only in 



VII 



1. The small A2 regime is in fact a subtle one since 



sec 
d = 

it is very sensitive to rare events and the connectivity 
properties of the Brownian path. 

We can furthermore deduce the asymptotic form of 
Fi(Ai) for Ai — > 00 and Ai — > by simple arguments as 
will be shown next. The behavior of F 2 (A 2 ) for large A2 
then follows from the fact that the large A2 behavior of 
eq. ( p2| ) should match the small Ai behavior of eq. (19), 
i. e. 1 - S- x F 2 {t/t%) ~ Fx(VC) for 1 < t£? < t « C . 



A. Limit Ai — > 00 

The particle leaves the site reached after t w typically 
in a time tt^ 1 . In order to explore the behavior for large 
Ai = t/t$ 3> 1 we may assume t t w 3> 1 and ask, 
which events give rise to a non-vanishing IT(t, S). These 
are rare events, where an unusual large inverse hopping 
rate T c g > t + t w ~ t has been encountered before the 
time t w has passed, or, said differently, before the particle 
has visited S(t w ) distinct sites. The probability that r e ff 
on one site is smaller than t is V(t) — J cLt c s p e s{T e s) ~ 
1 — t~ 6 /^~ a ' and the probability that at least for one of 
S = S(t w ) sites r c ff is larger than t is 1 — V s . Hence, 

fi(t, S) ~ 1 - V s = 1 - (1 - Q) s ~ Sit*,) Q(t) (23) 

where Q(t) = 1-V(t) ~ ^A 1 -"). Inserting S(t w ) from 
eq. (14a) we thus obtain 



F 1 (A 1 )~Ar 5 , S = — 

1 — ' 



(24) 



the probability ip(t)dt for the particle to leave the site 
reached after t w in a time interval (t, t + dt) can be esti- 
mated by 



!1!;lx > \ / dry 7~] J n 1 2 '2 



dT 2 T^ e 8 t 



=: <p(t,T m3X ) 



(25) 



Normalizing ip(t) on its support r m " x < t < r ma " we 
obtain 



i>(t, T m ax) = — 



Tmax) 



(26) 



J dt'<f>(t',T ma *) 



Since 1 — IL(t w + t,t w ) is the probability of the particle 
to leave the site reached after t w within the time interval 
[0, t), we expect 

l-IL(t w +t,t w )~ f dt'i>(t',T ma x(t w )) (27) 
Jo 

for t . fc,,) 3> 1 and t/tf^ 1 <§; 1. Performing the integrals in 
eq. (£5|) gives 



4>(t, Tmax) 



1 — a a-e 

t~ 

a-9 



1 max 



- 1 



(28) 



Since i/ r max ~ Ai, we have, in the limit Ai — > 0, to 
distinguish between a > 9 and a < 9, yielding 



4>{t, T m ax) 



1-q 
a-0 
1 -a 
0-a 



' max ) 



(29) 



a < 



Limit Ai 







Inserting this into eqs. (g6| g7|) and using ( |l7j ) then gives 
l-Fi(Ai) ~Af, Ai^O (30a) 

with 



In the limit Ai — ► corresponding to 1 < t < tjj 1 . the 
typical situations govern the decay of H(t w + t,t w ). In 
these typical situations, energies E <C -E m in(^M>) or Boltz- 
mann factors r 3> T max (/: u ,) do not matter. Let us then 
consider a jump of the particle from a site with energy 
E\ [Boltzmann factor r x = exp(— Ei/6)] to a neighboring 
site with energy E 2 [Boltzmann factor t 2 = exp(— E 2 /6)], 
where < Ei,E% < E m i n . On average such a jump oc- 
curs at a time t 2 /rl~ a . The probability for the particle 
to be at a site with energy in an interval (Ei, Ei + dEi) 
is ip(Ei) cxp(-Ei/9)dEi cx r-f 1_6 Vi<fri ("equilibrated 
initial site" ) and the probability for the energy of a 
neighboring site to be in an interval (E 2 ,E 2 + dE 2 ) is 
ip(E 2 )dE 2 oc t 2 x dr 2 ("random target site"). Hence 



a 

1 - 9 
1-a 



a > 



, a < 



(30b) 



As discussed above, by matching this small Ai behav- 
ior of Fi(Ai) to the large A2 behavior of F 2 (A 2 ) we fur- 
thermore obtain for a > 9 



F 2 (A 2 )~A^. 



A 2 — > 00 . 



(31) 



VI. TEST OF THE PEC 
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FIG. 5: First scaling function for a parameter set with a > 6: (a) H(t w + t,t w ) and (b) 1 — Tl(t w + t,t m ) as a function of 
t/t!^ 1 = Ai for parameters (9, a) — (1/4,3/8) and different dimensions d. The straight lines indicate the asymptotic behavior 
according to eqs. ( pi| , fjo| ). In order to make the graphs distinguishable, H(t w +t,t w ) and 1 — Tl(t w +t,t w ) have been multiplied 
by factors 4, 16 and 64 for d = 10, 100 and 1000, respectively. In plot (b) also t/tt^ 1 has been multiplied by factors 32, 4, 1 
and 0.5 for d — 1, 10, 100 and 1000, respectively. For d = 10 and d = 100 the symbols refer to the same waiting times as in 
Fig. §. For d = 1 the symbols refer to t w = 6 x 10 6 (+), 2 x 10 7 (o), 4 x 10 7 (♦), 1 x 10 8 (□), 3 x 10 s (■), 6 x 10 s (o), 2 x 10 9 (•), 
4 x 10 9 (a), and 1 x 10 10 (A). For d = 1000 the symbols refer to t w = 6 x 10 6 (♦), 1 x 10 7 (□), 3 x 10 7 (■), 8 x 10 7 (o), 2 x 10 8 (•), 
and 4 x 10 s (a). 



A. Scaling properties 

In order to test the PEC, we performed simulations 
of the model in d — 1,10,100 and 1000 by means of 
a continuous-time Monte-Carlo algorithm (see appendix 
[a| for details of the simulation procedure) . Averages were 
taken over 10 5 (for d = 1, 10, 100) and 10 4 (for d = 1000) 
energy landscapes. 

Figures [|a and||b show, respectively, H(t + t W: t w ) and 
1 — II (t + t w , t w ) as a function of the scaled variable 
Ai = t/t% (cf. eqs. (jl§||) for {9, a) = (1/4,3/8), i.e. 
a case where a > 9. As expected from the PEC and the 
scaling arguments outlined in the previous section, the 
data collapse onto master curves -Fi(Ai) for all dimen- 
sions d. In particular we find Fi(Ai) ~ A^ for Ai — > 
and i*i(Ai) ~ in agreement with eqs. (|30|a,b) and 
(p4|), respectively. 

Correspondingly scaled data for (9, a) = (1/3, 1/4), i.e. 
a case where a < 9, are shown in fig. ^. Again there is a 
good data collapse and the scaling functions -Fi(Ai) have 
the expected asymptotic behavior for small and large Ai . 
In p artic ular we now find -Fi(Ai) ~ A^ with e> from 
eq. ( [30a] ). 

Moreover, for a > 9, PEC and the scaling arguments 
predict the presence of a second time scale <2 = tf^ [cf. 
eq. (20)] and an associated gener alize d scaling t^, [1 — H(i+ 
t w ,t^y\ ~ F 2 {t/t^) [cf. eqs. (]2^Jl4a|)], The occurrence of 
this generalized scaling is verified in fig. |^ for d = 1, 
10 and 100. The master curves scale -£2^2) ~ A^ Q 



may note that the simulated data do not collapse at large 
A2, the deviations from scaling setting in at larger A2 for 
larger t w . The reason for this is that the limits t w — > 
00 and A2 = t/t%? — > 00 must not be commuted. One 
first has to take the limit t w — > 00 to obtain the scaling 
function ^(A2) and then has to consider the asymptotic 
behavior for large A2. 

The behavior for A2 — > is not so clear. For d > 1, 
the exact treatment of the PEC formula ( |l3"| ) yields 
-^2^2) ~ A^" also (with a smaller prefactor in the scal- 
ing law), while for d = 1, i<2(A2) ~ A2 seems to be 
correct. In fact, for d = 1 the data in fig. ^ are in fair 
agreement with a linear behavior. For d — 10, 100 the 
data also show some linear dependence for intermediate 
A2 values (see the dashed lines in fig. [?]), and in fact, 
such intermediate regime is predicted to occur based on 
the PEC formula (althou gh it s hould be less pronounced, 
see the discussion in sec. VII F ). For very small A2, how- 



for A 2 — > 00 as predicted by eq. (plj) . The critical reader 



ever, we expect i<2(A2) ~ A^ , although this asymptotic 
regime can not be identified in fig. |7| At best we can say 
that there is a change in curvature at very small A2. The 
true asymptotics, however, could not be obtained within 
reasonable CPU time. 

Accounting for the existence of two different time 
scales is very important to properly rescale the numeri- 
cal results. Had one assumed a single time scale t^, one 
would have obtained an approximate data collapse with 
an effective value of /i intermediate between /zi and fi2- 
This remark might be of importance for analyzing ex- 
perimental data: the assumption of a single time scale 
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FIG. 6: First scaling function for a parameter set with a < 6: (a) H(t w + t,t w ) and (b) 1 — U(t w + t,t m ) as a function of 
t/tw 1 = Aj f° r (0, a) — (1/3, 1/4) and different dimensions d. The straight lines indicate the asymptotic behavior according to 
eqs. ( p^ , po[ ). In order to make the graphs distinguishable II (tw + t,t w ) and 1 — H(t w +t,t w ) have been multiplied by factors 
4 and 16 for d — 10 and 100, respectively. In plot (b) also t/tt^ 1 has been multiplied by factors 32 and 4 for d = 1 and 10, 
respectively. The symbols refer to the same waiting time as for d — 1 in Fig. 0. 



could lead to a systematic underestimation of the true, 
asymptotic value of fj, (see the discussion in Q). 



B. Full comparison with simulations 

While the predictions of PEC regarding the scaling 
properties and the asymptotic behavior of the scaling 
functions are fulfilled, H(t + t w ,t w ) and its PEC equiva- 
lent TL(t + t w ,t w ) are different in quantitative terms. To 
see this, we have computed TV(t + t w ,t w ) numerically ac- 
cording to eq. ( [l3|) and compared it with H(t + t w ,t w ) 
obtained from the Monte-Carlo simulations. As shown 
in fig. ||the scaling functions associated with H(t+t w ,t w ) 
and H(t + t w ,t w ) differ by a factor in the scaling variable 
and the precise form in a transient region between the 
asymptotic regimes for small and large arguments. 

In summary we can conclude that all predictions of 
the PEC concerning the scaling properties can be cor- 
roborated by the simulations. The PEC thus turns out 
to be a powerful tool to uncover the mechanisms of aging 
in quenched random energy landscapes. 



VII. EXACT EVALUATION OF THE PARTIAL 
EQUILIBRIUM FORMULA 

In this part we show how the scaling arguments pre- 
sented in sec. ^ can be validated by an exact evaluation 
of the PEC formula ( |l3| ) in the limit of large S (large t w ). 
The reader who is not interested in these more mathe- 
matical derivations, may skip this section and proceed 
with the summary in sec. VIII. 



When we replace the denominator of eq. (13) by 

/ °° dX exp(— A Ylk=i Tfc ) an d notice that all but three 
random are uncorrelated with the Tj , r„ . appearing 
in the numerator of (13), we find 



IL(t,S) = S I dX (e~ XT f 3 g(t;X) , (32) 



where g(t; X) is defined by 
g(t;X) = 




-Ot _|_ —OL 



' I r l-c* 



t M 2 ^ 



(33) 



and 
/(*) 



9 dr 



a V a 



(34a) 



Here T(a,z) = J°° dte~ t t a ' 1 denotes the incomplete 
Gamma function. Note that < f(z) < 1 and that 
for z — * 



/(*) = i-r i 



z + 0(z' 



(34b) 



where T(a) = T(a, 0) is the Gamma function. For z — > oo 
(|argz| < 3tt/2) 



f(z) 

a z 



(34c) 
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In the limit S — > oo ( t w — > oo), the asymptotic form of 
eq. © is [see eq. @] 



Tl(t,S)~S d\e~ x s g(t;X) , (35) 
Jo 



where 



s — k s , k — r(i 



A. First scaling regime: Fi(Ai) 



(36) 



By substituting u — 5 1 / e A in eq. ( |35| ) and v = t/S 1 / 8 
in eq. (133) we obtain 



fi(t, S) ~ - / dae"" 



dv 



i r- 



-un/S 1 /" 



l r. 



1+0 



x exp — 1 6 



a , a 
1 ' T 2 



ts- 



(37) 

2(d-l) 



When keeping 



(38) 



fixed, (^3) nas a weu defined limit for 5 — ► oo. Before 
taking this limit, however, let us note at this point that 




< 2 = A 2 



FIG. 7: Scaling function i*2(A2) for the same parameter set 
as in Fig. |(| S[l — H(t w + t,t w )] as a function of t/tf^ 2 = A2 
for (0, a) = (1/4,3/8) and d = 1, 10, 100. The straight lines 
for l arg e A2 indicate the asymptotic behavior according to 
eq. (pll), while the solid (d = 1) and dashed lines (d — 10, 100) 
for small A2 correspond to a linear behavior (see text). The 
symbols refer to the same waiting times as in Fig. ^, with the 
additional waiting time t w = 3 x 10 6 . In order to make the 
graphs distinguishable, S[l — H(t w +t,t m )] has been multiplied 
by factors 10 and 5 for d = 1 and 10, respectively. 



if we would have considered all r nj in eq. ( |l3| ) to belong 
to the Brownian path, we had obtained 



fi(t, S) - 
2,1 



DO 

du e" 



dv 



1 T 4 



5-1/8 V" 



(39) 



instead of ©. Now, t>- e exp(-u e - w) E[j=i ^i" 
is an integrable majorant for the integrand in both 
eqs. ( |37| , ^9| ). Hence, by Lebesgue's theorem we obtain 
from both equations the same scaling function 



f\(Ai) = - / due~ u ~ j '^e~ uv 



f 



Ai 



2d 



(40) 



We can conclude that it makes no difference here whether 
we consider the path of distinct visited sites to have a 
one-dimensional or a compact topology or anything in 
between. As required by normalization, i*i(0) = 1. 



B. Limit Ai — > 00 

When transforming to variables £ = A^ 1 and 
C = Aj a 'v in eq. jic|), we find 

F l(Al ) = A- S ^ 6 - n d Ze-^ ni - a) t 6 x 
« Jo 

x / dtC e e- iC f 2d (C (1 ~ a) ) ■ 
Jo 

Since the integrand in the limit Ai — > 00 is an integrable 
majorant for all Ai, we can take the Ai — * 00 limit under 
the integral to obtain 



Fi(Ai) ~ Coo Ai 



-9/(l-a) 



for Ai — > 00 , (41a) 



where 



= 6 f°°_(K_ fad( C -(i- a )\ 



(41b) 



Note that this integral is well defined because of the 
asymptotic behavior of /(.) given in eqs. j34|b,c). 



C. Limit Ai -> 

Since iq(Ai) — » 1 for Ai — > 0, it is convenient to con- 
sider 



l-^i(Ai) 



du e 



l-f 



°dv 



2,1 



(42) 



10 



na w +t,t w ) 1 - n(t w +t,t w ) 




FIG. 8: Tl(t w +t, t w ) and l — H(t w +t, t w ) for (d, 6, a) = (10, 1/4. 3/8) as a function of the scaling variable t/t 1 ^ 1 from Monte-Carlo 
simulations (MC) and the Partial Equilibrium formula, eq. (y_3J) (PEC). The symbols refer to waiting times t w = 2 x 10 11 (+), 
6 x 10 11 (o), 3 x 10 12 (♦), 1 x 10 13 (□), 4 x 10 13 (■), 2 x 10 14 (o), 6 x 10 14 (•), 3 x 10 15 (a), and 1 x 10 16 (a). 



For 9 < a it follows from eqs. (|34|a-c) 1 — f(x) 2d = 
2dT{l-9/a)x e / a ip < (x), where (/?<(x) is a bounded func- 
tion, <^<(x) < M < for < x < oo, with y<(x) — > 1 for 
x — > 0. Hence, 



1 - ^(Ai) = 2dT 1 



- A 6 /" / due-' 1 
a/ « Jo 

dv 



— uv 

e (p< 



Ai 



1-a 



(43) 



Since M<i> e / a exp(— to — u e ) is an integrable majorant 
of the integrand, we find 

1 - Fi(Ai) - c< k\ /a , Ai -> 0, 6 < a , (44a) 
where 

2dOT(l-l) r°°. ... /- 00 cfo 



r(i 



due " , 

JO v 1 



^r(i-f) 2 r(i) 
r(i-e) 



(44b) 



For 6 > a the i; integral in eq. ( 43 ) would become diver- 
gent when taking the Ai — > limit in <£<(Ai/u( 1_a! )). We 
thus transform the v variable in eq. (|42]), w = Ai/v 1_Q , 
to obtain 



l-Fi(Ai) = 



(1 - a) « J 



du ( 



dw 



x [1 - f d (w)] exp 



mA 1 " 



(45) 



From eqs. (p^a-c) it follows for 8 > a that 1 - f 2d (w) — 
w(p > (w) , where (p > (w) is a bounded function, (p > (w) < 



My for < w < oo, with ¥>>(w) -> (0 - a)/2d0 for 
w — * 0. Hence, M> exp(— u 9 )k/ 1 ~ 6 ')/( 1_q ) is an inte- 
grable majorant of the integrand in ( p5| ) and we can take 
the Ai — ► limit under the integral, yielding 

l-Fi(Ax) ~c>Af^ Ai->-0, 6>a, (46a) 
where 



d.w 



u —[l-f d ( w )]. (46b) 



(i-a)r(i-fl)y w 1+ t 



D. Second scaling regime: J 7 2(A2) 

The second scaling regime is more difficult to extract 
from eqs. ( |32] , |33| ). We start by taking advantage of the 
normalization H(t — 0, S) ~ 1 in order to write 



1 - n(t, S) 




(47) 



This corresponds to eqs. ([37 39), if the number of neigh- 
boring sites belonging to the Brownian path is n, 2 < 
n < 2d (see the discussion in sec. IV). 

After the transformation v — » w = Ai $Z 7 - =1 r Q /ii'^ a ' 
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this can be rewritten as 

i-fl(t, s) eAl 

2d 



(1 - q)k J 



due— 



n 



2d 

j=i 



2(/ 



t2 -3=i r J dw (1 - e~ 



x exp 



We now decompose this expression into the sum of 



two parts corresponding to a decomposition Sf=i r j Q = 

2j=i r j Q + Sj= ra +i r j* m tne second line. We can use 
the symmetry with respect to {t^-} in order to replace 
Z)j=i r j" by nrf in the first part and $^ =n+1 rj* by 
(2d — n)T2 d in the second part. After the transforma- 
tion £ = «S -1 / e Ti in the first part, and the transformation 
C = S~ 1 / e T2d in the second part one can, for 9 < a, again 
use Lebesgue's theorem to perform the limit S — > oo for 
fixed 



A 2 = S a/e A! 



t 



S (l-2a)/0 

This yields 

S[l-II(t,S)] ~F 2 (A 2 ) 



1 

<a<- 



(48) 



(49a) 



F 2 (A 2 ) = nF 2 (1) (A 2 ) + (2d - n) F 2 (2) (A 2 ) , (49b) 



where 



fiK 1 2a )/ e r^i tjj 1 2a ^ e should increase with increasing 
t w . The function h(x) has the asymptotic behavior 

( (l-a)r(a-fl) , x->0+ 
~ < . (53) 

[ (l-a)r(l-6')a;- 1 , x -> oo 

The importance of the condition 9 < a now becomes 
clear, since for the £ integrand to be integrable in 
eqs. (BOlEa) for £ — > 0, a — has to be positive. 



E. Limit A2 — > 00 

After the transformation £ — > u = A 2 u^ 1_Q ^C Q , eq. ((5 
gives 



if>(A 2 



2 /-oo 



2 A 



a(l - a) k 7 u 1 - 9 /" 



dv 



(54) 



The limit A 2 — > 00 can then be taken under the integral, 
yielding [cf. eq. flU)] 



F 2 (1) (A 2 )~F 2 (2) (A 2 ). 
Using eqs. ( [l9b| , ^2] ) we finally obtain 



(55) 



p ca ^ J°°)A e / Q (00) _ 2dr(i - |) r(^) 

^2(A 2 )~c 2 A 2 , c 2 = _ . (5b) 

As required by matching, c 2 °°^ = c< [cf. eq. ( |44b| )] . 



F. Limit A 2 -> 



F 2 (1) (A 2 ) 



(1 — a) k Jo b" 



c 1 -^ 



MA^-^C) 



(50) 



fc(ss) ee x^" 1 r dw ^-^ W ) e -(-/«)^-/ (51) 



w 



and 



^ 2 (2) (A 2 ) 



(1 — a) k Jo u 



N/, ( «j- u ' /■ 00 dC/i(A 2 u 1 - Q C Q ) 




g 2 A 2 /a r°° due~ ue r°° dvh{v) 
a(l - a) k Jo u 1 - 9 /" Jo w / Q 



r(i-g rrf^) 
r(i 



■AS 



(52) 



The restriction to a < 1/2 in ( (48| ) follows from the fact 
that in order for the scaling regime to be relevant, t 2 ~ 



Since h(x) is monotonously decreasing with x for x > 
0, we can replace h(A 2 u 1 - a C a ) by h(0) = (l-a)T(a-9) 
in eq. (|50|). Hence, 



where 



c 2 



F 2 (1) (A 2 )~c 2 A 2 



] T[a-6) f°°due 



(57) 



dCe" 



■10 u Jo s v ' 

r(a-g)r m- 2 ("^) s 



1 - a 



r(i 



(58) 



This means that for small A 2 , F 2 (A 2 ) should be domi- 



nated by F 2 (A 2 ), except for d > 1 (where F 2 (A 2 ) = 
2-F 2 ^(A 2 )). For a given (mean) value of n we thus find 



F 2 (A 2 )^4 0) A^, 4°) s (i "U<»>. (59) 



This A 2 — * limit of the generalized scaling form depends 
on the number of neighbors being considered to belong 



12 



to the Brownian path (see the discussion in sec. IV). 
According to eq. (p7|), there can occur an intermediate 
regime, where ^(A^) depends linearly on A2. This in- 
termediate regime seems to be more pronounced for the 
"true dynamics" (see Fig. ^) than for the dynamics pre- 
dicted by the PEC formula @. 



VIII. SUMMARY AND CONCLUSIONS 

We have studied aging within the framework of a sim- 
ple hopping model mimicking a system that performs 
thermally activated transitions between the deep free en- 
ergy minima of its configuration space. Based on general 
arguments from the statistical theory of extremes we have 
chosen the free energy density of states to exhibit an ex- 
ponential tail. In order to effectively quantify the influ- 
ence of the initial and target site on the energy barrier to 
be surmounted during a transition, we introduced a pa- 
rameter a, < a < 1, in the hopping rates that turned 
out to strongly influence the aging properties. We have 
found that generically, subaging occurs in these models, 
an effect related to the multiple visits of deep traps. We 
have also found that different time scales, corresponding 
to different scaling regions, appear in these models. 

These aging properties can be understood from a par- 
tial equilibrium concept (PEC) that, despite not being an 
exact quantitative description, provides a powerful tool 
to study the scaling properties of the aging dynamics. 
Based on the PEC we first motivated the occurrence of 
subaging behavior and generalized scaling forms in terms 
of simple scaling arguments. We then presented a de- 
tailed analysis of the PEC formul a (|l3| ) and calculated 
the aging functions following from (|13[) and their asymp- 
totics exactly. With respect to the scaling properties the 
predictions could be confirmed by Monte-Carlo simula- 
tions in d = 1, 10, 100 and 1000 dimensions. 

The fact that even for d = 1000 the "quenched model" 
has aging properties different from the "annealed model" 
studied earlier in [j35| is rather surprising, since the num- 
ber of distinct visited sites S in d > 2 scales as the 
number of all transitions N between minima for large 
N, S ~ N. From this one tends to conclude that the 
system effectively explores a new minimum in each tran- 
sition, which would correspond to the annealed situa- 
tion. However, in the quenched situation one can imagine 
that there is always some local equilibrium established at 
the site with minimal energy reached after time t w , and 
this local equilibration effect slows down the diffusion 
in configuration space on all time scales, i.e. instead of 
S ~ N ~ i*' with 9' = min[l, 9/(1 - a)} in the annealed 
situation we have S ~ N ~ t 6 w in the quenched situation 
for d > 2, < 6 < 1. 

The existence of a local equilibrium around the "dom- 
inant" site with minimal energy after time t w does not 
imply that there must be a true equilibration on all vis- 
ited sites as it is assumed in the PEC formula (|l^). In 
fact, by studying the disorder averaged participation ra- 



tios Y q (t) = (EjPjM 9 ), wh ere Pj(t) is the probability 
for the system to be at minimum j at time t, we find that 
the PEC never becomes exact in the limit t — > 00, not 
even in d = 1 where each trap is visited an infinite num- 
ber of times. This behavior offers the possibility to define 
an effective temperature in the non-equilibrium aging 
regime, which enters a modified fluctuation-dissipation 
theorem, similar as it was found for mean-field spin glass 
models [ p8| . The role of an effective temperature in the 
landscape model considered here will be discussed else- 
where p0| (for recent progress in our understanding of 
this problem coming from MD simulations, see ]5l|). 

With respect to the applicability of the analysis out- 
lined above the question arises, whether the characteris- 
tics of the aging dynamics can be worked out also for gen- 
eral hopping rates not exhibiting the specific form given 
in eq. (||). For the PEC to be applicable, the system 
should have the tendency to approach equilibrium (that 
truly exists only for 9 > 1), so that one may require 
the jump rates to obey detailed balance. It is then in- 
deed straightforward, by using the simple arguments pre- 
sented in sec. |y|, to extract all characteristic time scales 
tj(t w ) and to predict the scaling properties. A necessary 
ingredient for this procedure to work correctly, however, 
is th e rob ustness of the scaling relation S(t w ) ~ [cf. 
cqs. (14a,b)]. Preliminary results indicate that for "phys- 
ical choices" of the jump rates (meaning that the depen- 
dence of the energies of the initial and tar get s ite on the 
saddle point energy is reasonable), eqs. (14a,b) always 
hold true. 

Moreover, it is possible also to consider some random 
distribution of the form of the jump rates (as it is ex- 
pected to occur when the dynamics in configuration space 
is mapped onto a jump process by means of some quan- 
titative analysis), and to work out the aging features of 
such more realistic models. It turns out then, that in 
principle infinitely many aging regimes can exists in the 
two-time plane t,t w > 0. A thorough discussion of these 
issues, however, is beyond the scope of the present work. 
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APPENDIX A: TECHNICAL DETAILS OF THE 
MONTE-CARLO SIMULATIONS 

We use the standard continuous-time Monte Carlo al- 
gorithm as discussed in detail e.g. in |i"9f to simulate the 



stochastic process defined in sec. [I A. A special problem 



arises for large dimensions d > 1, where it is not possi- 
ble to save the energies within a hypercube of even small 
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linear dimension. To resolve this problem, we use hash- 
maps, as, for example, the hash_map template provided 
by the Standard Template Library of ISO-C++. 

The hash-function should be computable quickly and 
at the same time the sites x being encountered must be 
mapped to different hash-values as often as possible. For 
dimensions d > 10 we found 



fix) 



2d 

E 



(Al) 



This yields eqs. (14a,b) for a = 



We note that in the annealed model eq. (B2) remains 

valid, while r max ~ A 1 / 9 , leading to S(t w ) ~ t^^ 2 for 
1 < d < 2, S(t w ) ~ t e w /\nt w for d = 2, and S(t w ) ~ t° w 
for d > 2. Due to our discussion in sec. Ill, one can 
replace 9 by 9' = min( 1,0/(1 — a)) in these formulae to 
obtain the behavior for < a < 1. 



2. One dimension (0 < a < 1) 



to do a good job. For large d the RAM consumption of 
the computer programs is the limiting factor when trying 
to access longer times t w . While simulations for d < 10 
can easily be performed on workstations, for d > 100 
computers with 4GB of RAM and more are necessary. 



APPENDIX B: CONNECTION BETWEEN THE 
"NUMBER OF DISTINCT VISITED SITES" AND 
THE WAITING TIME 

For a = 0, i.e. the trap model, we use a scaling argu- 
ment discussed in [f|7| to derive the behavior of S(t w ) for 
large t w . Then we show by a finite size scaling argument 
that in d = 1 the behavior of S(t w ) for a = is not 
expected to change for < a < 1. Furthermo re, w e give 
general arguments for the invariance of eq. (14q) with 
respect to a for all d. Finally we confirm eqs. ( |14a| ,b) by 
Monte-Carlo simulations. 



1. Trap model (a = 0) 



Let us consider a finite chain with L sites and site 
energies distributed according to eq. ([!]) . Then the mean 
square displacement (Ax 2 (t)) of a particle performing a 
random walk on this chain with the hopping rates given 
in (||) is expected to scale as 



{Ax 2 (t)) 



D(L)t for {Ax 2 (t))^>L 2 



for large L. The diffusion coefficient D(L) can be written 
as 



D(L) 



L' 2 



£.7 = 1 (P. 



(eq) 



(B5) 



w 



3,3+1, 



with p 



(eq) 



= exp(-/?^)/ELi ex P(-/3£fc) 



Tj/J2k T k- The denominator then reads, using (g), 
Y%=inY,j=i(TjTj+i)- a . Since (( Tj r j+1 )- a ) exists (for 
a > —9), the second sum gives a contribution oc L for 
large L, while the first sum has no finite average and is 
dominated by the maximum T max ~ L 1 / 9 . Thus we find 



After N 3> 1 transitions of the system, the typical 
elapsed time t w is 



JY 



i—1 fx 



N J drrp(T)^N[r max (S(N))}^ , (Bl) 



where g^N) is the typical number of Tj falling in some 
interval r M - Ar M /2 < 77, < 17, + At^/2, S(N) is the 
typical number of distinct visited sites after N jumps 
and r max (S) is the typical maximal r obtained after the 
system encountered S distinct sites. 
Since r max (S') ~ S 1 ' 6 and (see e.g. |1) 



N d ' 2 , 1 < d < 2 
S(N) - <( N/lnN , d = 2 
N d>2 



(B2) 



we find 



s [d+(2-d)8]/de j 1 < d < 2 
/„■(>') - ^ STi/^lnS d = 2 . (B3) 

SVfl , d > 2 



D(L) 



-1-1/ 



(B6) 



At the crossover time t^, where (Ax 2 (t x )} changes its 
behavior in eq. (B4), we obtain from continuity 



tl - D(L)t x ~ L 2 . (B7) 

This implies t x ~ L 2 /^ and ^ - L 1+1 / e yielding 

29 



V = 



1 



(B8) 



Since S(t w ) - (Ax 2 ^)) 1 / 2 , we finally obtain for < a < 
1 



(B9) 



in agreement with eqs. (14a,b) 



3. General arguments 

The very physical difference between the trap model 
(a = 0) and models with weighted rates (a > 0) is the oc- 
currence of forward-backward jump correlations. When 
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the system jumps from a site with low energy to a site 
with energy close to zero (such energies are most likely), 
it has high tendency to jump back for a > 0. More gen- 
erally, when the system enters a region of connected low 
energy sites, it will, before escaping this region, perform 
more and more jumps between the low energy sites the 
larger the value of a is. Once it left the region, it again 
has high tendency to jump back to it. 

One may regard a cluster of sites with deep energies 
and the surrounding shell of sites with higher energy as 
a "supertrap" . On a coarse-grained level with respect to 
time, the particle performs "superhops" between these 
supertraps. Then the essential difference between the 
a = and the a > situation disappears, since there are 
no increased backward jump correlations between the su- 
perhops. One thus suspects that a only rescales the time 
t w of the S(t w ) relation but does not change its expo- 
nent. Indeed this is what we have shown more explicitly 
in d — 1 in the previous subsection, and there is no reason 
why the argument should fail in higher dimensions. 

It is worth to note that the same arguments also apply 
to II(i + t w , t w ), if one generalizes it to a quantity II n (£ + 
t W it w ) which is defined as the averaged probability that 
the system after a waiting time t w does not leave a region 
of radius n in configuration space. Clearly, U(t+t w ,t w ) = 
IIo(i + t w ,t w ), but for larger n only superhops should 
lead to a decrease in II„(t + t w ,t w ). In fact we found 
that for n > 1 II n (t -I- t w ,t w ) shows normal aging, i.e. 

II»(t + t w ,t w ) ~ F[ n \t/t w ) for a > 0. 
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FIG. 10: Number of distinct visited sites S(t w ) for (9, a) 
(1/6, 1/4) and different dimensions d. 



(6, a) = (1/4,3/8) and (1/6,1/4), respectively. In all 
cases S(t w ) shows the behavior predicted by eqs. (14a, b). 



APPENDIX C: U(t, S) IN THE LIMIT S -> oo 

To derive the large S limit of eq. (pi 



Monte— Carlo results 



U(t,S)=S I dX (e- Ar ) S 3 5 (i;A) , 



o 



(CI) 
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we need to consider the small A limit of 
ddr 



-At 



-1+0 



= l-T(l-9)X e + - -A + 0(A 2 ) 

1 — C7 



where </?(A) is a continuous function with tp(X) — > T(l — 0) 
for A — > 0. Furthermore, (e _Ar ) has an upper bound 

e- aX " for A > 0, 



<e--) < 



— a A 



(C2) 



with some constant a, < a < 1, being independent of 
A. To proof this for A > 1 we compare 



FIG. 9: Number of distinct visited sites S(t w ) for (6, a) 
(1/4,3/8) and different dimensions d. 



—a A 



with 



= / dr 



i\ e T - i+d e- a ^y 



Figs, y and |l0| show S(t w ) from Monte Carlo simu- 
lations in d — 1, 10, 100, and 1000 for parameter sets 



/oo 
dr [Br- 1 -* e-*] 
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For 9 < a < 1 and r > 1 the first integrand is larger than 
the second integrand for all A > 1, thus the first integral 
is larger than the second integral. Since e~ aX for fixed A 
is strictly monotonously decreasing with a this rem ains 
valid also for < a < 9. To proof property ( |C2| ) for 
A < 1 we note that for small argument e it is 



0(e) = l-T(l-9)e e 

<l-ae e = e- a ^ +0(e W ) 



since r(l — 0) > 1 for 9 < 1. We can deduce that it 
exists a finite interval (0, Ao] where eq. (|C2|) holds. Be- 
cause (e~ AT ^ is strictly monotonously decreasing it is 
(e~ A ° T ) < 1. When choosing a = -ln(e~ A ° r ) > it 
holds that 



<e"^> < 



-an A 



for < A < 1 



and the proof of eq. (C2) is complete 



With the transformation A — » u = S 1 / e \ eq. (CI) gives 



rr(/.,s'i = I du cxp | i") u 



s 



When using eq. (|C3) we can estimate 



, M > (-*,(„«/,!?) 



S-3 , 
< exp I —a — - — u 



< exp 



k 

[6 

[7 
[8 

[9 

[10 
[11 

[12 

[13 

[14 

[15 



for S > 6 and from eq. (|33|), 

S^V'gfaS-V'u) < S 1 - 1 ' 6 



6dT 



-S- L,a UT 



< 01'(1 



[16 
[17- 
[18 
[19 



Thetie estimations will allow us Lo use Lebesuue's U10- 



orer 
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1 when cuiisideiing the liiniL 3 —y uu in eqs. (pTlyfl 



Thus we can write asymptotically 

/•oo 

Jo 



[21 

[22 
[23 



S I d\e- x kS 
la 



»(*;A) 



(C3) 



where the shortcut k = T(l — 9) [cf. eq. (J36|)] has been 
used. 
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